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Abstract 

The renormalization-scale (fx) dependence of the two-loop inclusive semileptonic b — > cl~vi decay rate is 
shown to be significant in the pole mass scheme, and the decay rate is shown to be poorly convergent in 
the MS scheme. Three-loop contributions to the decay rate are estimated by developing Pade approximant 
techniques particularly suited to perturbative calculations in the pole mass scheme. An optimized Pade estimate 
\^ , of the three-loop contributions is obtained by comparison of the Pade estimates with the three-loop terms 

determined by renormalization-group invariance. The resulting three-loop estimate in the pole-mass scheme 
^\ ' exhibits minimal sensitivity to the renormalization scale near fi = 1.0 GeV, leading to an estimated decay rate 

of 1927r 3 r(6 —¥ cl~Vf)l (G 2 F |V" C (,| 2 ) = 992±217GeV 5 inclusive of theoretical uncertainties and non-perturbative 
effects. 

o ; 

^CL,; 1 Introduction 

| The CKM matrix element \V c b\, which parametrises one of the sides of the unitarity triangle, can be extracted from 
the inclusive semileptonic decay rate V (B — > X c £~De). From a theoretical perspective, the inclusive process has 
the advantage that non-perturbative contributions are controllable; hence, an accurate perturbative determination 
of the b — > c£~vg decay rate is of value in obtaining \V c b\ from data. 

Complete two-loop calculations of semileptonic b — > c decays exist at the end points of the lepton invariant 
mass spectrum (maximal and zero recoil) [jjj and at an intermediate kinematic value Jg]. From these explicit 
calculations, the total semileptonic decay rate at two-loop order has been estimated [|2|7 In this present work, 
we extend these results to generate an estimate of the three- loop contributions to the b — > cl~vg decay rate via 
renormalization-group and Pade-approximant methods. 

In Section || we demonstrate that the pole mass scheme used in [2] has better perturbative behaviour than the 
MS scheme for the b — > c semileptonic decay rate, a result which is somewhat surprising since the MS scheme is 
better behaved within calculations of the b — ► u semileptonic decay rate || . The renormalization scale dependence 
of the two-loop b — ► c rate in the pole mass scheme is extracted using renormalization-group (RG) invariance, and 
the strong scale dependence that is found serves to motivate an estimate of next-order (three-loop) effects. 

The scale dependence of the decay rate is an important component of the procedure for estimating three-loop 
corrections to the b — > c semileptonic rate. Using information obtained from RG invariance, Pade approximation 
methods are developed in Section || that are appropriate for estimating next-order terms within pole-mass pertur- 
bative calculations. An optimized Pade estimate of the three-loop constant coefficient (i.e. the non-logarithmic 
term) is obtained by finding the best agreement between Pade estimates and true values of the RG-accessible 
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three-loop coefficients of logarithms. The RG-determinations of these latter coefficients in conjunction with the op- 
timized Pade estimate of the constant coefficient together constitute a scale-sensitive estimate of the full three-loop 
contribution to the perturbative rate. 

This Pade estimate of aggregate three- loop effects allows the renormalization-scale dependence of the b — ► ct~vt 
decay rate to be studied. In Section a region of minimal scale sensitivity [Q is found in the resulting decay 
rate. This minimal-sensitivity scale is found to be close to the fastest-apparent-convergence renormalization scale 
at which the three-loop contributions vanish entirely ^| . The proximity of these two scales supports the validity of 
these scales for obtaining estimates of the three- loop perturbative b — > cl~T>i decay rate. Theoretical uncertainties 
in this estimate are considered in Section 5l 



2 Scale Dependence of the Two-Loop Rate 



In ref. g], the inclusive semileptonic b — > cl i>i rate is estimated to two-loop order in terms of renormalization- 
group- (RG-) invariant pole masses m b and m c . The rate T can be expressed in the form 



r = 



Glml \ Vc 



cb 



1927T 3 

with the RG-invariant form-factor 

F(r) = 1 - 8r - 12r 2 log(r) + 8r 3 



^Miogf ^ 



\m b m c 



(1) 



(2) 



preceding the perturbative series S whose two-loop contribution at /i 2 = m b m c has been estimated in |^] by 
combining explicit results at the end and intermediate points Jl], ^) of the decay spectrum:[j] 



= 1 - 1.67 



a s (^/m b m c 



.9 ±0.3) 



(3) 



If we assume central values a s (m T ) — 0.33 |j, Q and m& = 4.9 GeV Q, and if wc follow ref. §] in assuming that 
m c = 0.3m b , we find that a s evolves from /i = m T through the four-loop, four-flavour ^ (rif = 4) /3-function [^) 
to a s (^m b m c ) jix = 0.087, in which case S [0.087,0] = 1 - 0.145 - (0.067 ± 0.002). It is evident from (§) that 
truncation after two-loop order introduces theoretical uncertainty of order 8.5%(= 0.067/0.79). 

There are, of course, other sources of theoretical uncertainty. Ambiguities concerning the definition of a pole 
mass in the presence of nonperturbative (confinement) effects have led the authors of ref. Q to reparametrise the 
rate ([!]) in terms of "low-scale" masses obtained through additional phenomenology — we will ultimately compare 
our results to the reparametrised rate in Section |5[ In light of more recent work || ||, |To[ |ll[ specifying accurate 
pole mass values by relating pole masses to MS and T-scheme masses with three-loop precision, we take a more 
empirical approach to the utilisation of pole masses within inclusive semileptonic rates. 

Renormalization-scale dependence provides an additional source of theoretical uncertainty to any rate calculated 
via the series (||). If S varies for different choices of /i, a value for T extracted from any particular choice of fi 
(e.g., fj, = m b m c ) is compromised. The optimal value of /i has been argued to be the choice for which S has 



minimal sensitivity to the renormalization scale 



the point at which dS/dfj, — 0. For (H) to lead to a 



reliable estimate of the true semileptonic b — > c rate, one would not necessarily need to establish that /i = m b m c 
is such a "Principle of Minimal Sensitivity" (PMS) point, but rather that the rate calculated at this value of /i 
differs only inconsequentially from the rate calculated at the true PMS value of /i. Any discrepancy between rates 
calculated at = m b m c and at upms is a direct measure of the former rate's theoretical uncertainty arising from 
renormalization-scale ambiguities. 

1 The O (oig) coefficient quoted in (^) [A. Czarnecki, personal communication] has been slightly corrected from the value —8.4 ±0.4 
appearing in ref. [§] . 

2 The use of four contributing flavours is necessary since o s is referenced to four flavours in (tj) M. 
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The two-loop order renormalization-scale dependence implicit in the perturbative series within ([!]) may be 
parametrised as follows: 

S [x(fjt), L(/i)] = 1 + (o + aiL) x + (b Q + b 1 L + b 2 L 2 ) x 2 , (4) 

a Q = -1.67 , b Q = -8.9 ± 0.3 , x{n) = , L(/i) = log ( — — \ . (5) 



m h m c 

The renormalization-scale invariance of the all orders rate implies that 

" 2 ;£ = » < 6 » 

hence that, 

dS[x,L] tn o n % n a \ dS\x,L] 
0= dL - (f3 x 2 + (3 lX 3 + f3 2 x 4 + . . .) ( ? ) 

where the normalization of QCD /3-function coefficients is explicitly defined by 

fJ , 2 ^ = -((3 x 2 +Pix 3 +p 2 x* + ...) (8) 

with n f = 4 values (3 = 25/12, 0i = 77/24, and j3 2 = 21943/3456 (in the MS scheme). It is easily seen from (0) 
that the logarithmic coefficients within (Q) are 

at = b 2 = , bi = a [3 = -3.479 . (9) 

In Figure |[ we have plotted the mass- and scale-dependent portion of the rate (Q), 

as a function of //. The Fig. [I] curve is obtained by assuming rtib = 4.9 GeV ||, m c — 0.3ma Q, and the central 
value bo = —8.9 from the estimate Q given in (^). Substantial scale dependence is evident from the figure: the 
rate increases monotonically with ji 1 flattening out somewhat for larger values. Moreover the curve exhibits no 
PMS point (i.e., extremum), and at \i = ^/m c m,\, = 2.68 GeV yields a rate 10% smaller than the (still increasing) 
rate at fj, = 8 GeV, the flattest portion of the curve shown in the figure. Thus, the two-loop pole-mass calculation 
of the b — > cl~ vi rate exhibits at best only poorly-controllable dependence on the choice of renormalization scale. 
Such scale dependence may compromise any subsequent "low-scale" mass expression devolving from the two-loop 
/i = ^/mfcTO c pole-mass rate. 

Similar scale dependence and even worse apparent non-convergence characterise the two-loop order pole-mass 
calculation of the b — > u£~&£ rate, thereby motivating a recasting of the calculation in terms of the MS running 
6-quark mass ||. For four contributing flavours, the two-loop relationship between pole and MS running quark 
masses is ficfl 



m pole = m(ji) 



I + * G4)) + ( io - 3919 + f * G4) + 1 * (J® ) ) i2 « 

in) 



By substituting this relation into ([!]) , (^) , and Q , we find that the fully MS version of the mass- and scale-dependent 
portion ([l0|) of the b — > c&~vt decay rate (using the central value 60 = —8.9) is 



r ^=^)^(^|) {l + (5.00 + 6Iog(^)]. < -„O 



x (12) 

J "' :! - 1 71 1 ' sfcj) + 4 " bs (4m) + f ^ G&O ] l2,rt ■ 
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Note that F (rn 2 /m 2 ^ remains RG invariant if m c — > Tn c (fi) and to;, — > mj(fi). If /i 2 = mb{^)m c {^) and if we 
continue to assume that m c (/i) = 0.3mb(/i), then 

-— = m 5 b (^)F (0.09) [1 - 1.02 x (//) + 18.2.x 2 (//)] , fj, = V(L3m b (fJ,) (13) 

The convergence of this MS perturbative series is even more ill-behaved than its pole-mass version (^). Moreover, 
the scale dependence of dl2| ) is shown in Figure || to be even more pronounced than that of the same rate in the 
pole-mass scheme (Figure]!]) — Figure |^ displays a rate which decreases with fi with no apparent PMS point. 

Thus the MS approach, which substantially improves the perturbative series within the semileptonic b — » u rate 
|J, fails to improve the pole-mass expressions (|l|-||) for the semileptonic b — > c rate. If this latter rate is to be 
utilised to extract an estimate of \V c b\ from the inclusive B — ► X c i~vi branching ratio, there is evident value in 
having an estimate of next-order corrections in order to obtain some control over renormalization-scale dependence. 
This three-loop order contribution to S [x(/i), L(/j,)] is necessarily of the form 

AS 3L [x(ji), L(jjl)] = [co + ciL(ji) + c 2 L 2 (fi) + c 3 L 3 ( M )] x 3 {^) (14) 

The RG equation (fa) implies that 



=aix + (bi - a /3 ) x 2 + (2b 2 - /3 ai) x 2 L + [a - 2/3 b Q - a p{\ x 3 + [2c 2 - 2f3 b 1 - foai] x 3 L 
+ [3c 3 - 2f3 b 2 ] x 3 L 2 + O (x 4 ) . 



(15) 



The set of results (||) are evident from the requirement that O(x), O (a; 2 ) and O (x 2 L) terms in ([TH]) separately 
vanish. The vanishing of subsequent terms in (15) implies that 

d = 2b /3 + Oo/3i , c 2 = a o 0l , c 3 = . (16) 

For ao and &o as given in (|s|) and three massless flavours (as appropriate to phenomenology devolving from a s (m T ) 

i), 

Cl = -42.4 ± 1.3 , c 2 = -7.25 , c 3 = . (17) 

The coefficient cq, however, is RG- inaccessible to these orders of perturbation theory, and requires a direct three- 
loop calculation. In the absence of such a calculation, we estimate cq in the section which follows via asymptotic 
Pade approximant methodology, in much the same way as in a prior estimate ]l3[ ] of the three-loop contribution to 
the b — > u£~Dg decay rate. 



3 RG/Pade Estimate of cq 

Consider a perturbative field-theoretical series with N + M known terms: 

S(x) = 1 + R lX + R 2 x 2 + . . . + R N+M x N+M + . . . . (18) 

The set of known coefficients {Ri, R 2 , ■ ■ ■ , Rn+m} is sufficient to determine in full the N + M coefficients charac- 
terising an [JV|Af] Pade-approximant to the series S: 

1 + a x x + a 2 x 2 + . . . + a N x N 

S[N\M]{X) = ——r 2 — TT ■ 19) 

1 + bix + b 2 x l + . . . + b M x M 

The coefficients {a!,... , ajy,6i,... are obtained by the requirement that the power-series expansion of 

S[n\m]( x ) recovers the N + M known coefficients within (18), the series S(x). The next O (x JV+M+1 ) term in 
this power series is a Pade approximant prediction for the first unknown coefficient Rn + m+i- For example, if only 
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the next-to- leading order coefficient R\ is known, one can use this coefficient to construct a [0|1] approximant to 
the series (|l8|), 

S m (x) = = I + R lX + Rjx 2 + . . . , (20) 

1 — JX\X 

for which R\ is the predicted value of R2 : 

4 0|1] = R\ ■ (21) 

Similarly, if R\ and R2 are known (corresponding to two sublcading orders of perturbation theory) , one has enough 
information to construct a [1|1] approximant to the series (^8|): 

1 + f R 2 i-R2 \ x 

S m( x ) = r £ =l + R 1 x + R 2 x 2 + ^x 3 + ... . (22) 

1 - ^-x Ri 

The predicted value for the unknown coefficient R3 is 

4 1|1] - Rl/Ri (23) 

In general, one can always use an [JV|M] approximant ( |l9| ) to predict the first unknown series coefficient 
Rn+m+i- Such predictions have accuracy which increases as N and M increase. For perturbative field-theoretical 
series (characterised by asymptotic Rn ~ N\ behaviour) the accuracy of such predictions has been argued by Ellis, 
Karliner and Samuel to satisfy the relative error formula [[l4| 

adWM] p[N\M] _ ntrue i\t\aM 
^^N+M+l _ n N+M+l ^N+M+l - A 

^+m+i " <+m+i (iV + M + aM + 6) M [ ' 



where A, a, and b are constants to be determined. Of particular interest are the relative errors obtained from (24) 
for the predictions ( ^l| ) and (pjj) 

fl i- fl 2 - d (25) 

R2 l + (a + 6)' 1 5j 

iz^ = (26) 

i? 3 2 + (a + b) ■ 1 j 

Denoting the constant a + b = k, we can eliminate the other constant A within (^5|,^6|) and solve for R3 algebraically 
to obtain the improved estimate 

_ (2 + k)Rl 

U3 (l + k)Rl + R X R 2 ' ( ' 

The assumption fc = has been utilised in prior applications to predict successfully the third subleading contribution 
to the QCD MS 0- function 0, the DRED SQCD 0- function the MS /3- function for massive scalar field theory 
|lq ], as well as to obtain estimates of such contributions for a number of processes calculated in the MS scheme: 
SM @ and MSS M M Higgs -> 2 gluon decay rates, 6 -> ||, VT+iy- ->■ ZZ 0], and the QCD static 

potential function |l9|. 

The choice fc = 0, however, is ill-suited to pole-mass calculations. To see this, consider the series [l4|) within 
(Q) characterising the pole-mass expression for the b — > c£~ &g decay rate 

S[x, L] = 1 + (oo + aiL) x+ (b + b x L + b 2 L 2 ) x 2 + (c + c x L + c 2 L 2 + C3L 3 ) x 3 (28) 

x = , L = log 1 



m b m c 
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We have already seen that ao = —1.67, a± — 0, 6 = —8.9 ± 0.3, b% = aoPo — —3.48, b 2 — 0. The above series is in 
the form of (Eq) with Ri and R 2 respectively identified with ao + a\L and bo + b\L + b 2 L 2 . If As = 0, we see from 



(27) that 



_ 2 (M- g| £ , 26, 

ao(ao + ^o + &ii) ao ao 

Comparing this expression to the form ( fTi)) anticipated for the third subleading order of S[x,L], we necessarily 
obtain the following predictions for the RG accessible coefficients C2, c±: 

b 2 

c 2 = 2^ = 2a /3 2 , (30) 
ao 

d = 2 (26 -<%) — = 4/3 &o ~ 2/3 a 2 (31) 
a 

It is evident from ( |l6| , p"7| ) that these predictions are quite poor; (|3(]) is double the true value for c 2 , as obtained 
in ([l7|), and c\ is also badly overestimated [for values (||) and the central value bo = —8.9, eq. (|3l]) implies that 
Ci = —85.8, in contrast to the correct (RG) value c\ = —42.4]. 

For a given choice of k, estimates of the coefficients a characterising the third subleading order have been 
obtained for correlation functions [M and the b — > ui~Dg rate |l^] by moments of R3, as estimated in (p7|), over 
the entire ultraviolet region (e.g., / (mi,m c ) > 1 for the case at hand). These moments are then equated to 
corresponding moments of R3 = Co + c\L + c 2 L 2 + C3L 3 in order to obtain values for {co, ci, C2, C3}. If we define 
w = mbiric/iJ, 2 (L = — logw), the moments 

1 

Nj = (j + 2) J w 3+1 R 3 (w)dw (32) 


can be obtained using (|27j) for the integrand i?3, which becomes a function of u> for our pole mass case by virtue 
of the w-dependence of R 2 : R\ — ao, R 2 = bo — 61 log w. After (numerical) computation of the values of Nj, the 
coefficients Cj are obtained by equating such values to the corresponding integrals 

1 

Nj = (j + 2) J w 3+1 (cq - c\ log w + c 2 log 2 w - c 3 log 3 w) dw . (33) 


In particular, we see that 

7V_i =c + ci + 2c 2 + 6c3 (34) 
113 

No=c + -ci + -c 2 + -c 3 (35) 
12 2 

^Vi = c + -ci + -c 2 + -c 3 (36) 
^2 = c + Jci + ic 2 + ^c 3 . (37) 



If k ^ — 1, such a procedure is seen to lead to a non-zero value of C3, in contradiction to the result C3 = (16) 
necessarily following from application of the RG equation (0) within the pole mass scheme. For the case k = — 1, 
however, the result ( p7|) collapses to the naive estimate (|2^), which in the pole mass scheme is necessarily a degree-2 
polynomial in L: 

4 k -V = M = 1 ( bo + b lL f = b l L 2 + 2^L + h l (38) 
Ri ao ao ao ao 
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The moment procedure described above then reduces to fitting the degree-3 polynomial R3 = C3L 3 +c 2 L 2 +c\L+cq to 
the degree-2 polynomial in (|38|). Such a fit necessarily reduces to equating the powers of L in these two expressions. 
We thus find that C3 must be zero, and that 

b 2 

c 2 = — = ao/? 2 (39) 
a 

ci = 2^ = 2p bo (40) 
ao 

co = ^ (41) 
a Q 

Equation ([59]) is in exact agreement with the result ([l6|) obtained via RG-invariance from (0). Equation ( |40| ) 
represents the first (and dominant) contribution to the RG-determination of ci in ([l6]). For the central value 
bo = —8.9 (and n/ = 4) the ( |40| ) prediction ci = — 37. f is not far from the true value ci = —42.4. The accuracy 
of these results provide some support to the naive estimate cq = —47 .4 obtained via ([id]) for the RG-inaccessible 
coefficient cq. 

We can further improve our estimate for cq by finding the value of k within (^) which, when used within the 
integrand of (|||) to match the moments Nj to (Q-^) , most closely reproduces the true values of C2 and c\ , as 
determined by RG methods in ( [l6| ) and (0). In such a procedure we obtain estimated values for C3, C2, Ci, and cq 
that depend explicitly on k: 

1 27 32 

c (fc) - -gJVLi(fc) + 4iV (fc) - yiVi(fc) + yiV 2 (fc) (42) 

3 189 

ci(fc) = - 32AT (fc) + — JVi(fe) - 647V 2 (fc) (43) 

ca(fe) = -y^V-i(fc) + 38^) - ^iV x (fc) + ip7V 2 (fc) (44) 

2 32 
c 3 (fc) = giV_i(fc) - 8iV (fc) + 18JVi(A) - yiV 2 (fc) (45) 



where 

1 

iV i (fc) = (i + 2) / 







(2 + fc)(& -b 1 log W ) 3 
(1 + fc)ag + a (b - 61 \ogw) 



dw (46) 



We then use the explicit values for c\ and c 2 obtained from (JT^) by RG methods to optimize the sum of the squares 
of the relative error of estimated values of c\ and c 2 , 

A ^ _ f c 2 {k) -aof3j \ 2 | / ci(fc) - (2f3 b - a ffi) \ 2 ^ 



with respect to k. For the set of values ao = —1.67, bo = —8.9, (3q = 25/12, (3\ — 77/24, we find a clear minimum 
of A(/c) at fc = —0.94, as evident from Figure [|, consistent with the case made in the preceding paragraph for an 
optimal k value close to k = —1. Corresponding k = —0.94 values for the moments Nj are obtained numerically 
via (E6|) 



7V_i = -106.3 , N Q = -74.92 , N x = -66.17 , N 2 = -62.12 (48) 

and these values lead via (|4^)-(^)to the following estimates for the third subleading order coefficients: 

c 3 = 2.0 x 10~ 4 , c 2 = -7.68 , ci = -39.7 , c = -51.2 (49) 

These values reflect excellent agreement with the RG values (|l7|). The above estimate for the RG-inaccessible 
coefficient Co is only 20% larger in magnitude than the naive prediction ([ill), indicative of the internal consistency 
of the methodology. 



7 



CO 




- Cl 


-2c 2 - 


- 6c 3 = - 






1 


1 


3 


CO 


2 C1 


--c 2 


"4 C3 = 




= Nt- 


1 


2 


2 


CO 


3 C1 




"9 C3 = 
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CO 


4 C1 


~8 C2 


"32 C3 = 



We conclude by noting that four separate estimates of Co can be obtained from (|34j)-(37) by substituting into 
these equations the optimal Nj values ( f48| ) as well as the true values C3 = 0, c\ = —42.4 and C2 = —7.25 as 
determined by RG invariance (n7h in the previous section. Solving each equation separately for Co, we find that 

49.4 , (50) 

-50.1 , (51) 

-50.4 , (52) 

= -50.6 . (53) 

The estimates ( pfj| ) ( p3| ) are remarkably consistent with each other and only slightly smaller in magnitude than the 
estimate in ( [49| ) obtained without RG inputs for {c\, C2 ,03}. 

In estimating the third subleading order of Jl4|), the series within the rate ([!]), we choose ( pl| ) as the most 
central of our cq estimates, in conjunction with the explicit RG determinations of c\ and C2 (C3 = 0). This set of 
values is obtained, as noted earlier, for the central value estimate Q of &o : 

b Q = -8.9 : {c ,ci,c 2 } = {-50.1, -42.4, -7.25}. (54) 

Precisely the same procedure can be used to obtain corresponding estimates of {co, ci, C2} for the extremes of the 
bo = —8.9 ± 0.3 range obtained in Q (see footnote |l|): 

&o = -9.2 : {co,ci,c 2 } = {-53.6,-43.7,-7.25} (55) 
b Q = -8.6 : {cd,ci,C2}S {-47.5, -41.2, -7.25}. (56) 

These estimates are obtained in precisely the same way as those of ([54| ): ci and c 2 are identified with their RG 



values via (g_7J), and cq is determined via (51) with No(k) [see Eq. (46)] evaluated at a value of k which minimizes 
A(/c) [see Eq. (ftfl)]. The coefficient c 3 = for all cases, as evident in the previous section from RG-invariance. For 
the 60 = —9.2 case, the minimizing value k = —0.94 is the same as for bo = —8.9. For bo = —8.6, the minimum 
value of A(fe) is found to occur when k = —0.93. 



4 Scale Dependence of the Three-Loop Rate 

If 60 = —8.9, the three-loop b — ► c£~&i inclusive estimated rate is given by (|l|) with 

S [x((j,),L(jj,)] = 1 - 1.67a; + (-8.9 - 3.479L) x 2 + (-50.1 - 42.4X - 7.25L 2 ) x 3 . (57) 

The coefficients of x 3 are those of ( |54| ) , as estimated in the previous section. Figure |] displays a plot of the fi- 
(scale-) dependence of the "reduced" three-loop rate ( |l0| ) with S[x,L] given by (|57|). The pole masses mi, and 
m c are assumed to be mt, = 4.9 GeV and m c = 0.3m,b = 1.47 GeV consistent with values used in Figure 
^ clearly displays a much flatter /^-dependence than the two-loop rate plotted in Fig. 0. In addition to this 
diminished dependence on the renormalization scale fi, the rate plotted in Fig. |J also exhibits a distinct minimum 
at fi — 1.0 GeV. At this PMS (minimal-sensitivity) value of /i, the successive terms of the series S[x,L] exhibit 
reasonable convergence: 

S [x(1.0 GeV), £(1.0 GeV)] = 1 - 0.258 - 0.049 + 0.020 (58) 

Fpj\£5 _ r(l.OGeV) _ 2Q47Q e y5 _ (59) 
K K 

Eq. (|9|) corresponds to the minimal-sensitivity value for the rate, as discussed in Section |2[ 

Note that the small three-loop contribution to ( |5S| ) can be tuned to zero by making only a small change in the 
choice of the renormalization-scale parameter \i. The value of /1 at which the three-loop term vanishes (i.e., the 
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value of fi at which the Fig. [I] two-loop and Fig. ^| three-loop curves intersect) corresponds to the renormalization 
scale associated with the "fastest apparent convergence" (FAC) of the series S[x, L\. This occurs at /i = 1.18 GeV: 

S [x(1.18 GeV), i (1.18 GeV)] = 1 - 0.226 - 0.058 + (60) 

r^ = r(i.i8GeV) =1Q5lGeV5 

K K 

It is striking that the FAC 0] and PMS 0] criteria predict virtually identical rates; moreover, a similar equivalence 
of rates obtained via these same two criteria is found for the estimated three- loop contribution to the b — > u(.~i>i 
rate In both semileptonic processes, the PMS and FAC momentum scales are comparably small. The PMS 
and FAC scales (1.0 GeV and 1.18 GeV) for b — ► cl~vi are respectively 37% and 44% of the logarithm reference 
scale y/mbm c = 2.7 GeV. For b — > u£~Dg, H PM s — 1-78 GeV and fi FA n — 1-84 GeV, numbers which are 42% and 
44% respectively of the logarithm reference scale mj, (m^) = 4.2 GeV || |l3| |. 

Although the estimated three-loop rate plotted in Fig. ^ exhibits much less dependence on the renormalization 
scale \i than the two-loop rate of Fig. [l], we anticipate the existence of residual /i-dependence as a consequence of 
the truncation of the series d57j ) after three-loop order. A way to eliminate much of this residual scale dependence 
is to "undo the truncation" by choosing an appropriate Pade approximant to the series (|57j). For example, a [2|1] 
approximant to the series ( |57| ) that reproduces its power series to O (a; 3 ) is 

where 

1 ' 8.9 + 3.48L 1 ; 

Similarly, the known series terms ( |57j) are also reproduced in the power series of the [1|2] approximant 

^ /T , 50.1 + 42.4L + 7.25L 2 + 1.67 (8.9 + 3.48X) , . 

Ei(L) = (67) 

v ' 1.67 2 + 8.9 + 3.48L v ' 

(8.9 + 3.48L) 2 - 1.67 (50.1 + 42.4L + 7.25L 2 ) 

E2{L) = 1. 67 2 + 8.9 + 3.48L (68) 

Dx (L) = Ei (L) — 1.67 (69) 

In Figure ^| we have superimposed plots of the reduced rate (10) using 



1. Eqs. (|^) and (^|), the two-loop version of S[x, L] leading to the reduced rate also plotted in Fig. 

2. Eq. (|57j), the three-loop version of S[x, L] leading to the reduced rate plotted in Fig. ||, 

3. S I 2 ' 1 ', as determined via (^32|-|65|), and 

4. Si 1 ! 2 !, as determined via <MM). 



The vertical scale of Figure || is magnified compared to that of Figs. [I] and || in order to accentuate the differences 
between reduced rates obtained for each of the above scenarios. We observe from Fig. [|that both Pade-approximant 
versions of the rate coincide after \i = 1.5 GeV and are considerably flatter than the rate devolving from the three- 
loop version of S[x, L]. Indeed, the three-loop reduced rate is itself quite stable, increasing slowly from 1066 GeV 5 to 
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1180 GeV 5 as /j, increases from 1.5 GeV to 8 GeV. Nevertheless, the two Pade approximant versions of S[x, L] vary 
only minimally over the same range of /z, increasing from 1058 GeV 5 at /i = 1.5 GeV to 1079 GeV 5 at /J, = 8 GeV. 
Thus Pade-improvement of the three-loop rate virtually eliminates the residual scale dependence of the naively 
truncated expression. Such use of Pade approximants to eliminate residual scale dependence is also evident in Fig. 
3 of ref. |l3|] for the b — > u£~Di rate, and has been previously discussed in the context of the Bjorken sum- rule [^lj 
as well as in more general terms ^2|. Of particular interest, however, is the convergence of all four curves in Fig. 
p| to virtually the same PMS/FAC point. This convergence lends further support to the PMS/FAC estimates (p9|), 
(61) for the reduced rate. 



5 Discussion 



The entire analysis presented in Section [| can be repeated using the extreme values bo = —8.6 and bo = —9.2, as 
estimated in Q (again, see footnote |l|), utilising (|55|) and (|56|) for the appropriate determinations of three-loop 
coefficients in conjunction with the known values of ao and b\ (|9j). One finds the uncertainty in bo is reflected in a 
±14 GeV 5 spread in the PMS/FAC value 1050 GeV 5 for the reduced rate. 

Other sources of theoretical uncertainty arise from a s (m T ) = 0.33 ± 0.02 |, Q], mf le = (4.9 ± 0.1) GeV {§, 
and the error that may occur in estimating cq. We estimate that Pade determinations of cq are subject to errors 
comparable to those of Pade determinations of the RG-accessible coefficients c\ and c 2 ; e.g. 
7%. If we are conservative and estimate the uncertainty of Co to be double that of ci, 

Sc 



)/c{ 



RG I 



CO 



14% 



(70) 



the corresponding uncertainty in the reduced rate is ±38 GeV 5 . Consequently, our estimate of the purely pertur- 
bative three-loop order b — ► ct~vi rate, as defined by (10), is 



ppe 



K 



(1050 ± 14 ± 44 ± 115 ± 38) GeV 5 



(71) 



where the listed theoretical uncertainties respectively devolve from the uncertainty in bo, ct s (m T ), m^ ole , and cq. 



Nonperturbative (NP) contributions to the rate may be extracted from Eq. (5. 
following additional contributions to the series S[x, L]: 



AS 



NP 



Ai + 3A 2 
2mg 



6 1 



A 2 



where the form factor F is given by (g) and where 

-0.5 GeV 2 < Ai < , A 2 = 0.12 GeV 2 . 
This additional NP contribution entails a ~ 6% reduction in the reduced rate (|7l| 



r 

K 



K 



.5 ± 5.5) GeV J = (992 ± 217) GeV 



if) and in the intermediate step of (p 



of |23J , and correspond to the 



(72) 

(73) 
(74) 

have been combined 



where the independent sources of uncertainty in 
additively. 

If we identify the predicted b — > cl~i>i decay rate with the inclusive semileptonic process B — > X c l~vg \l — e or 
//, but not their sum], we can then relate the aggregate theoretical uncertainty in (|^) to the concomitant theoretical 
uncertainty in the determination of \V c b\'- 



\Vcb\ = 



192tt 3 SBR(S X c £-& e ) 



-,1/2 



GpT B [(992 ± 217) GeV 5 ] 



(75) 
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To factorize experimental and theoretical uncertainties, we employ recent central values for the average B lifetime 



t b [29] and the B — > X c e v e branching ratio [E6| to rewrite (75) in the following form: 



The first factor in (f76|) reflects the summed theoretical uncertainties in ((73), which are separately broken down in 



(71). We have been conservative in identifying and assessing the magnitude of each such independent source of 
error — the theoretical uncertainty estimated from a (truncated) two-loop calculation of the b — > cl~vt rate should 
be larger than that obtained by us in ( |76| ) . 

Note also that we can compare our T pert /K = 1050 GeV 5 central value estimate ( |7l| ) for the reduced rate 
(exclusive of NP effects) with the corresponding estimate one would obtain using low-scale masses 



T pert fm 2 c 

= rrCi F — | 

K h \mt 



1-1.14 K ' - (3. 5 ±0.3) V 



(77) 



If we utilise the low-scale mass values fh% = 4.64 GeV, m c = 1.25 GeV, as quoted in || from ref. [E4J and find 
via devolution from a s (m T ) — 0.33 that a s fy 'inline) /tt = 0.091, we observe that the rate predicted via ( |77j ) is 
1097 GeV 5 , an answer in surprisingly good agreement with our 1050 GeV 5 three-loop estimate in the pole-mass 
renormalization scheme. 

We conclude by noting that the RG equation (|L |?]) may be used to determine additional higher-order corrections 
to the decay rate ([!]). The leading- log corrections to all orders in x are determined by the one- loop /3-function; 
next-to-leading-log corrections to all orders in x are determined by the two-loop /3-function etc. Although we 
have made use of RG-invariance to O (x 3 ) in present work, it is in fact possible to incorporate these logarithmic 
corrections to all subsequent orders. The full exploitation of RG-invariance within perturbative-QCD expressions 
for physical processes is presently under study. 
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Figure 1: Renormalization scale (fi) dependence of the two-loop reduced rate T/K in the pole mass scheme. 
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Figure 2: Renormalization scale (/i) dependence of the two-loop reduced rate T/K in the MS scheme. In the MS 
scheme, mb(n) is obtained from the four-loop nf = 4 anomalous mass dimension [fl2"| using mj, (to?,) = 4.2 GeV [§) 
as a reference value. 
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Figure 3: The quantity A measuring the sum of the relative errors in the Pade estimate ( |4^j44| ) of the RG- 
accessible three- loop coefficients {ci, C2} plotted as a function of the error- formula parameter k (|27|). The location 
of the curve's minimum represents the value of k leading to an optimized Pade estimate. 
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Figure 4: Rcnormalization scale (//) dependence of the three- loop reduced rate Y/K in the pole mass scheme. The 
PMS point is represented by the local minimum of the curve. 
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Figure 5: Scale dependence of different estimates of the reduced rate in the pole scheme. The solid curve represents 
the three-loop estimate presented in Fig. and the dotted curve represents the two-loop estimate also presented 
in Fig. [l]. The [1|2] and [2|1] Pade approximants obtained from the three-loop estimated rate are represented by 
the dashed curves which overlap almost completely above /j, — 1.5 GeV. The PMS point is represented by the local 
minimum of the three-loop curve, and the FAC point occurs at the intersection of the two- and three-loop curves. 
Note the convergence of all the estimates near the FAC/PMS points. 
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